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The Higgs boson production pp (pp) — bh + X via bg — bh at hadron colliders, which may be an 
important channel for testing the bottom quark Yukawa coupling, is subject to large supersymmetric 
quantum corrections. In this work the one-loop SUSY-QCD corrections to this process are evaluated 
and are found to be quite sizable in some parameter space. We also study the behavior of the 
corrections in the limit of heavy SUSY masses and find the remnant effects of SUSY-QCD. These 
remnant effects, which are left over in the Higgs sector by the heavy sparticles, are found to be 
so sizable (for a light CP-odd Higgs and large tan 3) that they might be observable in the future 
experiment. The exploration of such remnant effects is important for probing SUSY, especially in 
case that the sparticles are too heavy (above TeV) to be directly discovered in future experiments. 


14.80.Cp, 13.85.Qk,12.60.Jv 


I. INTRODUCTION 


Searching for the Higgs boson is the most important task for the Fermilab Tevatron collider and the CERN Large 
Hadron Collider (LHC). Among various Higgs production mechanisms, those induced by the bottom quark Yukawa 
coupling are particularly important because in some extensions of the Standard Model (SM) such a coupling could 
be considerably enhanced and thus the production rates can be much larger than the SM predictions. The Minimal 
Supersymmetric Standard Model (MSSM) [1] serve as a good example of such extensions, where the coupling of 
the lightest C P-even Higgs boson (denoted by h) to the bottom quark is proportional to tan 8 [2] and thus can be 
significantly enhanced by large tan 8. 

In the production channels of the Higgs boson via its coupling to the bottom quark , the process pp (pp) > bh + X 
via bg — bh was recently emphasized in Ref. [3]. The advantage of this process over the production via bb — h [4], the 
dominant production channel via the bottom quark Yukawa coupling, is that the final bottom quark can be used to 
reduce backgrounds and to identify the Higgs boson production mechanism [5]. And compared with the production 
via gg, qd — hbb [6], the production rate of pp (pp) — bh + X is one order of magnitude larger. So the production 
pp (pp) — bh -- X may be a crucial channel for testing the bottom quark Yukawa coupling. 

If the MSSM is indeed chosen by Nature, then the prediction of the cross section for the production pp (pp) > bh-- X 
[7] must be renewed with the inclusion of SUSY quantum corrections because, like the process of the charged Higgs 
boson production pp (pp) —^ tH~ + X [8,9] and the relevant Higgs decays [10-12], the SUSY quantum corrections to 
this process may be quite large. In this work we study the one-loop SUSY-QCD corrections to this process, which is 
believed to be the dominant part in the SUSY corrections. 

It is well known that the low-energy observables in the MSSM will recover their corresponding SM predictions when 
all sparticles as well as M4 (the mass of the CP-odd neutral Higgs boson) take their heavy limits. If only some of the 
masses take their heavy limits, e.g., all sparticles are heavy but M4 is light, then large remnant effects of SUSY may be 
left over in the physical observables of the Higgs sector. The study of these remnant effects will serve as an important 
probe for those heavy SUSY particles [11]. Such kind of study will be performed for the process pp (pp) > bh + X in 
this work. After deriving the SUSY-QCD corrections to this process, we will examine the behavior of the corrections 
in the limit of heavy SUSY masses. When MA is light and the sparticles take their heavy limits, we find the large 
remnant effects left over by the SUSY-QCD in such a Higgs production process. 

This paper is organized as follows. In Section II we present our strategy for the calculation of the one-loop SUSY- 
QCD corrections to the process pp (pp) — bh + X. In Section III, we scan the parameter space of the MSSM to 
estimate the size of the SUSY-QCD corrections. In Section IV, we study the behaviors of these corrections in the limit 
of heavy SUSY masses. The conclusion is given in Section V and the detailed formula obtained in our calculations 
are presented in the Appendix. 


II. CALCULATIONS 


At high energy hadron colliders, the incoming b-quark is generated from gluons splitting into nearly collinear bb 
pairs. When one member of the pair initiates a hard-scattering subprocess, its partner tends to remain at low pr and to 
become part of the beam remnant. Hence the final state typically has no high-pr b-quarks. When the scale of the hard 
scattering is large compared with the b-quark mass, the b-quark is regarded as part of the proton sea [13]. However, 
unlike the light quark sea, the b-quark sea is perturbatively calculable. If the scale of the hard scattering is ji (for the 
scale we use ji to distinguish from the SUSY parameter p), the b-quark distribution function b(z, ji) is intrinsically 
of order os(íi) log(fi/mp). As ù approaches m, from above, b(z, jt) — 0; while as j becoming asymptotically large, 
os(ft) log(ft/my) approaches order of unity and one needs to sum terms of order o? (A) log" (A/m) into the b-quark 
distribution function to yield a well-behaved perturbation expansion in terms of a, [13]. In this case, the b-quark 
distribution function becomes of the same order as the light partons. The main uncertainty of b-quark distribution 
function comes from that of gluon distribution function which is about 1096 [14]. 

The subprocess gb — bh occurs through both s-channel and t-channel shown in Fig. 1(a,b). The spin- and color- 
averaged differential cross section at tree-level is given by 
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where 4 and f are the usual Mandelstam variables, œs(ñ) is the running strong coupling, and mp(j:) is the running 
bottom quark mass [4]. « represents the mixing angel between the two CP-even Higgs boson eigenstates and 8 is 
defined by tan 6 = v3/v, with v1.2 denoting the vacuum expectation values of the two Higgs doublets [2]. In Eq.(2.1), 
we use the MS running mass of the b-quark rather than the pole mass to take into account large QCD logarithm 
corrections to the vertex hbb [15]. The SM prediction of the cross section is recovered when setting | sin a/ cos 9| = 1 
[3]. Throughout the calculations we neglect the b-quark mass except in the b-quark Yukawa couplings. 

The one-loop Feynman diagrams of SUSY-QCD corrections are shown in Fig. 1(c-r). In our calculations we use 
dimensional regularization to control the ultraviolet divergences in the virtualloop corrections. For the renormalization 
of strong coupling constant gs, we employ the MS scheme [16]. As to the hbb Yukawa coupling, at one loop level to 
O(as) it is given by 


= CD SQCD 
Ühbb = Ghbb + gS sb: gi , (2.2) 


where gabb denotes the one-loop coupling, gn, is tree-level coupling, bas is the radiative correction from pure 
QCD [15], and bgp CP is the one-loop SUSY-QCD contribution [11]. In determining bapa + Ogres one needs 


= QCD S D 
the counter-term of the vertex hbb, whose general from is given by ghb»( — + m, ) with ôm, being the 
counter-term of the b-quark mass defined by m? = mp + ómy (m$ is the bare mass). ómy is determined by requiring 
my to be the pole of the one-loop corrected b-quark propagator [15,11,17]. One major difference between Carag and 


gg pct is that the former contains large logarithms a, lo -Æ of O(1) and thus one needs to introduce MS runnin 
hbb Sin 8 
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mass my(ft) to absorb leading logarithms a? log (E) [15]. An extensive discussion about this issue in the MSSM 
was provided in Ref. [17]. 
The one-loop SUSY-QCD contribution to the amplitude of gb — bh can be written as 
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where Pr. rR = (1 F 35)/2, T^ = A*/2 with A? being the Gell-Mann matrices, and k, p; and pə are the momentum of 
the incoming gluon, incoming b-quark and the outgoing b-quark, respectively. gs and my should be understood as the 
running ones in Eq.(2.1). The coefficients C; arise from the loops and are given explicitly in Appendix A. We have 
checked that all the ultraviolet divergences canceled as a result of renormalizability of the MSSM. 


di di 


FIG. 1. Feynman diagrams of gb — bh with one-loop SUSY-QCD corrections: (a,b) are tree level diagrams; (c — e) are 
one-loop vertex diagrams for s-channel; (f — k) are one-loop vertex diagrams for t-channel ; (l— n) are the box diagrams; (o— r) 
are self-energy diagrams. 


The differential cross section of gb — bh with one-loop SUSY-QCD corrections is given by 
dê _ dè? | d(AG) 
dí dt dÊ 


where the first term is the tree-level result given in Eq.(2.1) and the second term is the one-loop SUSY-QCD corrections 
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The cross section of gb — bh is then given by 
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where £44,, = 0 and £44, = —84- m? with mj; denoting the Higgs mass. In order to avoid collinear divergence in 


Eq.(2.6) and to enable the outgoing b-jet to be tagged by silicon vertex detector at the Tevatron and the LHC , we 
require the transverse momentum of the outgoing b-jet to be larger than 15 GeV and apply a rapidity cut || < 2.5 
for the LHC and |m] < 2.0 for the Tevatron. 

The total hadronic cross section for pp (pp) — bh + X can be obtained by folding the subprocess cross section & 
with the parton luminosity 


o(s) = / dr a G(s = sT), (2.7) 


where ro = m2 /s and s denotes the pp (pp) squared center-of-mass energy. dL/dr is the parton luminosity given by 
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where f? and f? are the b-quark and gluon distribution functions in a proton, respectively. In our numerical calcula- 
tions, we used the CTEQ5L parton distribution functions [18]. We did not distinguish the factorization scale Q and 
the renormalization scale (4, and assumed ñ = Q = mp. The scale dependence of our results will be briefly discussed 
in the proceeding section. 

The process pp (pp) — bh + X has been extensively studied [3] in the framework of the Standard Model. Its 
cross section is found to be at the order of 1 fb for the Tevatron and 100 fb for the LHC, and the next-leading-order 
(NLO) QCD correction can enhance the production rate by 50% ~ 60% for the Tevatron and 2096 ~ 40% for the 
LHC [3], depending on the applied cuts and the Higgs boson mass. We will incorporate such QCD corrections in our 
calculations for the production rate o/osm. 


III. NUMERICAL RESULTS 


In this section we will perform a scan over the SUSY parameter space to show the possible size of the SUSY-QCD 
corrections. Before performing numerical calculations, we take a look at the relevant parameters involved. For the SM 
parameters, we took my = 80.448 GeV, mz = 91.187 GeV, m, = 174.3 GeV, mi (mp) = 4.2 GeV [19], sin? Oy = 0.223 
and as(mz) = 0.118 . We used the one-loop QCD running o(fi) and ma (A). 

For the SUSY parameters, apart from gluino mass, the mass parameters of sbottoms are involved. The sbottom 
squared-mass matrix is [2] 
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mj = met mg mz (I5 — Qi sin? 6w) cos(28), (3.2) 
mi. = m% + m? + m2, Qy sin? Ow cos(28), (3.3) 
Xp — Ap E i. tan B. (3.4) 


Here ms and m are soft-breaking mass terms for left-handed squark doublet Q and right-handed down squark D, 


respectively. A, is the coefficient of the trilinear term H 1QD in soft-breaking terms and u the bilinear coupling of 
the two Higgs doublet in the superpotential. 72 = —1/2 and Q, = —1/3 are the isospin and electric charge of the 
b-quark, respectively. This mass square matrix can be diagonalized by a unitary rotation 


br, — cos 0p — sin 0j b 
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and consequently 6, and the masses of physical sbottoms b12 can be expressed as 


tan 20, = —3 E (3.6) 
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From Eqs.(2.1,2.5,2.6) we know that the cross section also depends on the Higgs mass, a and B, which can be 
determined at tree level by tan 8 and the CP-odd Higgs mass M4 [2]. Noticing the fact that both the mass and the 
mixing angle receive large radiative corrections when the SUSY scale is high above m; [20], we used the loop-corrected 
relations of Higgs masses and mixing angle [21,22] in the computation of cross section. In our calculation, we used 
the program SUBHPOLE2 [21], where two-loop leading-log effects of the MSSM are incorporated in the Higgs masses 
and the mixing angel, to generate Mp and a needed for our computation. The input parameters for this program are 
the mass parameters in the top sqaurk and sbottom sector, and M4, tan B and the heavier chargino mass my. 


We found that the usage of the loop-corrected relations of Higgs masses and mixing angle is indeed necessary. 
Comparing with the results obtained by using tree-level relations for Higgs masses and mixing angel, the size of 
SUSY-QCD corrections by using the loop-corrected relations is generally magnified from 30% to 20096!. We also 


checked that this conclusion is also valid for the SUSY-QCD correction to the Higgs partial width r(A — bb). 
Now we know the relevant parameters are 


mo, Mg, Mp, Åt, b, mg, ms, pi, Ma, tan B, (3.9) 


where Mz is the soft-breaking mass term for right-handed top-squark and A; the coefficient of the soft-breaking 


trilinear term H3QU. To show the main features of SUSY effects in pp (pp) — bh + X, we performed a scan over this 
ten-dimensional parameter space. In our scan we make no assumptions about the relations among these parameters 
to keep our result model-independent, but restrict the parameters with mass dimension to be less than 2 TeV. In 
addition, we consider the following experimental constraints: 


(1) u > 0 and tan f in the range 5 € tan 8 € 50, which seems to be favored by the muon g — 2 measurement [23]. 


(2) The LEP and CDF lower mass bounds on Higgs, gluino, stop, sbottom and chargino [24,25] 


mn = 114 GeV, m;, > 86.4 GeV, mj, > 75.0 GeV, mg > 190 GeV, mz?" > 67.7 GeV, (3.10) 
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where mie is the mass of the lighter chargino. 


30 | UN 
š ^  tanp- 8 
ex .* © tanp-30 
. + tanB=50 
"r Pu 4 
ek the 
NEA n 
E a; 
S ET ae 
m St Bt 
3 k 
m i +. 
ata 
3 ee t 
ost , 
4 
f. 
2 LER 
F 
+ Sip cote. 
ERI " 
L + 
* uL TIR TET] "m 
EUNT] 
1 r- T ig ad wd + 
0.9 " 
0.8 NI l l pal 


FIG. 2. The scatter plot of o? /a$ y versus Ma. 


It would be interesting to first scan over the allowed parameter space to figure out how large the production rate is 
enhanced in the MSSM. In Fig.2 we present the tree-level cross section relative to the SM prediction with the same 
Higgs mass. This ratio is independent of collider energy, but is dependent on the SUSY mass parameters since we use 
the loop-corrected relations of the Higgs masses and the mixing angle (see the second paragraph of Sect. II) . From 
Fig.2 one sees that the production rate in the MSSM can be significant larger than the SM prediction for a light M4; 
while for a heavy M4 the MSSM prediction approaches to the SM value. This character was first noticed in [26] and, 


'The main reason for such an enhancement is that for a large SUSY scale, the dominant term of the SUSY-QCD correction 
is proportional to cot a + tan (see for example, Eqs.(4.2, 4.4, 4.6, 4.9) in the proceeding section), whose value obtained by 


using the loop-corrected relations of the Higgs masses and the mixing angle is generally larger than that by using the tree-level 
relations. 


as a result of this character, distinguishing the lightest MSSM Higgs boson from the SM Higgs boson in the large 
Ma limit will be very difficult. When SUSY-QCD corrections are added, this character remains unchanged for heavy 
sbottoms(see following discussions). From Fig.2 one also finds that there exists the possibility (although very rare) 
that the MSSM cross section is suppressed to be below the SM value [27]. In this case, the SUSY-QCD corrections 
will play a more important role in Higgs phenomenology at colliders [28]. 

Now let us scan over the allowed parameter space to show the possible size of the SUSY-QCD corrections relative 
to the tree-level value. In our numerical evaluation, we found the relative correction is insensitive to collider energy. 
'The difference of the results between the LHC and the Tevatron is at the level of parts per mill. In Fig.3 we present 
the SUSY-QCD corrections to the Cross section. 
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FIG. 3. The scatter plot of the SUSY-QCD correction Ac /a? versus MA for the Tevatron. The difference of the results 
between the LHC and the Tevatron is at the level of parts per mill. 


Fig.3 manifests three features of SUSY-QCD corrections for large tan 8. The first one is that the correction size is 
enhanced by tan 9 and thus can be quite large. The second one is that for MA lighter than 500 GeV, the correction 
tends to be negative. The third one is that for large M4, the correction may be positive and the maximum value 
seems to be independent of the value of M4. These features can be explained as follows. 

For the correction size larger than 2%, the dominant contribution of the SUSY-QCD correction is from the loop 
corrections to the vertex hbb and the contribution of the box diagrams is much smaller for the parameters satisfying 
the constraints in the paragraph following Eq.(3.9) ?. As a result, the correction behaves like (which is similar to the 
SUSY-QCD correction to the vertex hbb in case of heavy sbottoms [11]) 


Ac N Mw M$w 


(3.11) 


where Mew and M; denote the electroweak scale and the typical mass of sbottoms, respectively. Cı and C2 are 
functions of m; , m; , mg, p and Ap, but independent of MA. It is found that in general C4 is negative and C% is 
positive and either C1 or both C, and C5 are enhanced by large tan B. For a light M4 compared with mg, the first 
term of the RHS in Eq.(3.11) is dominant and hence the cross section tends to be negative. While for a large MA , 
the second term is dominant and the cross section tends to be positive. So the behavior of Eq.(3.11) can explain the 
features of Fig.3. 

From Fig.3 we noticed that in some corners of parameter space the one-loop SUSY-QCD contributions to the cross 
section are comparable or even larger than the tree-level result and consequently, one must consider higher order 


corrections. In such cases, it is important to sum over the terms o7( M = )* to all orders of perturbation theory by 


using an effective Lagrangian approach [17,29]. 


?In the large limit of SUSY mass parameters discussed in the proceeding section, we have checked that, even for the correction 
size far smaller than 1%, the dominant contribution still comes from the corrections to the vertex hbb. 
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FIG. 4. The SUSY-QCD correction óc [o9 and the cross section c/o si versus the mass of the produced Higgs boson. Solid 
curves are for the LHC and the dashed for the Tevatron. For Ac/o?, each solid curve overlaps with the corresponding dashed 
one due to the tiny difference. 


Next we study the dependence of the SUSY-QCD correction óc /c? and the cross section normalized by the SM 
prediction, o/agm, on the mass of the produced Higgs boson, which can be directly compared with experiment 
results and hence is much informative. In such a study we assumed a common value (Msysy) for all input SUSY 
mass parameters and, considering the fact that mp is insensitive to M4 for M4 > 150GeV [22], we fixed the value of 
Ma. Then through varying the value of Msysy, we obtain the different mass value of the Higgs boson. 

The dependence on mp is illustrated in Fig.4 for tan 8 = 50. (Note that mj, can vary only in a small range since 
it is stringently upper bounded in the MSSM.) In this figure and also in the following figures showing o/osy, we 
also incorporated the conventional QCD corrections [3] into both the MSSM and the SM cross sections. As pointed 
out earlier, the SUSY-QCD correction 0c /o? is not sensitive to collider energy. The difference of the results between 
the LHC and the Tevatron is too small to be visible, as shown in the upper part of Fig.4. But for o/ogy the 
difference between the LHC and the Tevatron is visible since the QCD corrections are significantly different for these 
two colliders [3]. 

We also studied the dependence of the production rate on the renormalization scale ji. (As pointed out earlier, we 
assume that the factorization scale is equal to the renormalization scale.) We found that such a scale dependence is 
significant in some parameter space. For example, for the LHC with MA = 300 GeV and mp = 120 GeV, the ratio 
c(ü)/c(mj5) is 0.93 for & = ma/2 and 1.03 for i = 2mp. Such an uncertainty is comparable with the uncertainty 
from the b-quark Yukawa coupling (rny(rmy) = 4.2 + 0.2) and the partron distribution function (about 10%). 


IV. BEHAVIOURS OF SUSY-QCD CORRECTIONS IN DECOUPLING LIMITS 


To study the behaviors of the SUSY-QCD correction in the large limit of SUSY mass parameters, we consider 
four typical cases as in Ref. [11] where the decoupling property of SUSY-QCD correction to the coupling of Abb is 


analyzed. To qualitatively understand the feature of each case, we present the approximate formula in the limits, but 
in practical numerical calculations we use the complete one-loop expressions. 

About the inputs of the SUSY parameters, there are several differences between our work and Ref. [11]. The first 
one is that in [11] the tree-level relations for the Higgs masses and the mixing angle were used, but in our calculations 
we use the loop-corrected relations. As discussed earlier, using the loop-corrected relations leads to a significantly 
different correction. The second one is that in our analysis we considered the experimental bounds in Eq.(3.10). This 
will rule out some parameter space which have been considered in Ref. [11]. The third one is that in cases B and D 
we also require A+» to be large since large A; is favored by the Higgs mass bound. 
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FIG. 5. The SUSY-QCD correction Ac/o? versus Ms in Case A. For each fixed value of M4, the solid curve (for the LHC) 
overlaps with the corresponding dashed one (for the Tevatron) due to the tiny difference. 


(1) Case A: All SUSY mass parameters except M4 are of the same size (collectively denoted by Ms) and tend to 
heavy, i.e., 


mo ~ mi e mg e Ay c» A, ~ mg ~ po Ms. (4.1) 
In this case the SUSY-QCD correction behaves like 
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where the first term in the RHS corresponds to the first term in Eq.(3.11)? and the rest corresponds to the second 
term in Eq.(3.11). The striking feature of this case is that for very large Ms, the correction approaches a nonzero 


3In the MSSM the tree-level relation for Higgs masses and mixing angle predicts the following relation: cota = — tan 8 — 


2 4 
ae tan B cos 23 + OE) and at loop level, the gap between cot a and — tan @ is generally enlarged. 


constant, and this remnant effect of SUSY-QCD corrections is enhanced by large tan 8. This feature is illustrated in 
Fig.5. From Fig.5 one also finds that the SUSY-QCD correction in this case is negative and sizable * for a light MA 
and a large tan 6. 

In Fig.6 we show the loop corrected cross section normalized by the SM prediction. From this figure we see that for 
MA of several hundred GeV, although the tree-level cross section in the MSSM can be reduced by large SUSY-QCD 
corrections, an enhancement of several times over the SM prediction can still be expected due to the fact that the 
tree-level b-quark Yukawa coupling in the MSSM is significantly enhanced by large tan 8 for light Ma. This large 
enhancement shows a very weak dependence on Ms. So we can conclude that up to the next-leading order, a light 
Ma is still able to make the MSSM cross section larger than the SM prediction. 
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FIG. 6. o/osm versus MA in Case A for the LHC (solid) and the Tevatron (dashed). 


(2) Case B: Mg Mg, Mp and Az,» (collectively denoted as Ms) is much larger than u, mg and MA (collectively 
denoted as M ), i.e., 
Mö p ~ Ate ~ Ms > ms c p o Mar M. (4.3) 
In this case the SUSY-QCD correction behaves as 
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“Tf the correction is too sizable (say exceed 50%), higher order corrections are also important and need a proper treatment 
[17,29]. 


From this expression we see that in the large Mg limit, the SUSY-QCD corrections decouple rapidly as M?/M 3 and 
the decoupling behavior is slowed down by large tan 8. The characters of this case are shown in Fig.7. So we see that 
even with a fixed light M4, the process still does not have remnant SUSY-QCD effects if the gluino mass and p are 
also kept light. 
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FIG. 7. The SUSY-QCD correction Aoc /o° and o/asm versus Ms in Case B. The solid curves are for the LHC and the 
dashed for the Tevatron. For Ac / o°, each solid curve overlaps with the corresponding dashed one due to the tiny difference. 
(3) Case C: Only the gluino mass gets much larger than other SUSY parameters (collectively denoted as Mg) °: 
mg > Mò ü,Ď ~ Ara ~ ur Ma ~ Ms. (4.5) 


In this case the SUSY-QCD correction behaves as 


Ac 2a, |2Ms M2. Mscota m? 
Ara t ii- ee 
a sia) — ugs MÀ 
Ms tan B mz cos Bsin(a+ 3) mf tan? Bcota (4.6) 
M; M sina M3Ms ) 


2 
'The main character of this case is that as gluino mass gets large, the correction drops very slowly like + log x 
g S 


which was also observed in Refs. [11,12]. Again, like other cases, the size of the correction is enhanced by large tan f. 
In Fig.8 we show the dependence of the SUSY-QCD correction and a/agy on the gluino mass. Note that in this 


In this case the Higgs mass bound requires Ms to be much larger than electroweak scale. 
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case we found that tan 8 = 8 cannot satisfy the experimental bounds in Eq.(3.10) for Ms = 600 GeV. In Fig.8 the 
correction size is significantly smaller than those in Case A and B. The reason is that here a large M4 is chosen so 
that tan 6 + cosa is suppressed (see footnote 3). 

Let us explain the origin of es slowness of the decoupling i in this case. Such slowness of the decoupling arises from 


the first term in Eq.(4.6), i.e., D <=“ (tan 8 + cot a) log = (note that the Ms in the factor E is u and the one in the 


VEI 
logarithm is squark mass Mj). As M; gets much larger than Mj and p, i decreases but " r increases. For the 
g a 


example shown in Fig.8, i.e., Mg is oe from 1 TeV to 5 TeV with fixed Mg = u = 600 GeV, the factor KA is 


decreased by 1/5 but the factor log Ma af? is increased by 4.16. Thus the slowness of the decoupling as gluino gets heavy 
is caused by the enlarged mass splitting between gluino and squark. Of course, since tan @ + cot o is proportional to 
2m —# tan cos 23 (see footnote 3), the contribution of the first term in Eq.(4.6) will be decoupled rapidly if MA gets 
Den 


0r : : 
L Case C J 
2003 - Ma-Mg-M =A, ,=u=M, = 600 GeV d 
(o 906 - tan p = 30 
» |" 
^l  .009 


t tanp = 50 
-0.12 


Tos 


` `1000 2000 3000 4000 5000 


Gluino Mass (GeV) 
FIG. 8. The SUSY-QCD correction Ac /o? and c/osw versus the gluino mass in Case C. The solid curves are for the LHC 


and the dashed for the Tevatron. For Ac/ o°, each solid curve overlaps with the corresponding dashed one due to the tiny 
difference. 


(4) Case D: One of the sbottoms and A; become heavy while other mass parameters (denoted as M) are fixed. 
We choose 


Mp ~ Mg ~ Ate > Mo ms u~ Ma~ M Mew (4.7) 
or equally 


m, ~ mig ~ Atb > mg n mj c Mgr Ur M > Mew, (4.8) 
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where mi. and mz, 12 are the masses of bottom-squarks and top-squarks, respectively. 
In this case the SUSY- QCD correction behaves as 


D|- (4.9) 


Cirig 
m2 sin a 3°W 


Ac 20, E 
b2 mj, 


M 2 i 
—— rw (tan B + cot a) 1+ logy „ ma cos p sin (a + 8) 
c 30 


2 mi 


'The main feature of this case is the correction decouples like TA log “+ — and this decoupling behavior is slowed down 
Mia 

by large tan B. Fig.9 shows the dependence of the correction and the cross section on mj (~ mi. ). Comparing with 

the results in case B, we see that the correction size decrease more slowly as mp gets heavy. 


0 
-0.03 — 
6 E 
~ -006 f 
o6 L 
<] L 
-0.09 — 
ARE QUE MANI 
400. ^. ` abo 1200 1600 2000 
148 — 
8 a 
© 346 [- 
< L 
o L 
134 — 
L Í ii | 4 L | L Í L 
400 800 1200 1600 2000 


FIG. 9. The SUSY-QCD correction Ac /o? and o/osm versus Mp in Case D for the LHC (solid curves) and the Tevatron 
(dashed curves). For Ac/c?, each solid curve overlaps with the corresponding dashed one due to the tiny difference. 


From the above analyses we see that when MA is fixed and all other SUSY mass parameters get large, the SUSY- 
QCD lefts over some remnant effects in the Higgs production process pp — bh + X at the hadron colliders. Note that 
for the remnant effects to be left over, y and gluino mass must be comparable with or larger than the masses of the 
sbottoms. The fundamental reason for such a behavior is that the couplings like hb;b; are proportional to SUSY mass 
parameters. 

We conclude this section by making a few remarks. Firstly, in our analysis we assumed mg, u > 0 and, as a result, 
all four cases have negative values of the correction. In the anomaly-mediated SUSY breaking scenario [30], a negative 
mg is predicted and in this case, the sign of the correction may be reversed. Secondly, it should be noted that in 
the calculations of such Higgs processes it is necessary to use the loop-corrected relations of the Higgs masses and 
the mixing angle since such relations can significantly affect the results. Thirdly, in case of a light My, although the 
SUSY-QCD corrections tend to reduce the cross section severely, the MSSM cross section can still be several times 
larger than the SM prediction. If the cross section of this process is measured in the future and found to be several 
times larger than the SM prediction, a light M4 is favored. 
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V. CONCLUSION 


In this work, we studied the one-loop SUSY-QCD quantum effects in the Higgs production pp (pp) — bh + X at the 
Tevatron and the LHC in the framework of the MSSM. We found that for a light M4 and large tan 8, the corrections 
can be quite sizable and cannot be neglected. We performed a detailed analysis on the behaviors of the corrections 
in the limits of heavy SUSY masses and found that when M4 is fixed and all other SUSY mass parameters get large, 
the SUSY-QCD lefts over some remnant effects in the Higgs production process pp — bh + X. Such remnant effects 
can be so sizable for a light M4 that they might be observable in the future experiment. The exploration of such 
remnant effects is important for probing SUSY, especially in case that the sparticles are too heavy (above TeV) to be 
directly discovered in future experiments. 
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APPENDIX A: EXPRESSIONS OF FORM FACTORS 


Before presenting the explicit form of Cis, we define the following abbreviations: 


$— (pi + k)’, t= (k — po)’, (A1) 
d ir 1 : 
123 = Jav Oy F cos 6b), bi 2 = Er O» + sin 6b), (A2) 
AT ES (ar Fs br)’, AT = (ar + br)’, Ar = aj us bi, (I = 1,2), (A3) 
AT, = aja; + bby +arbj + bray, AT, = aja; + bib; — arb; — bray, (A4) 
Br, = ajaj — brbj — arbj + braz, BP, = arag — bib; + ajo — bray, (A5) 
1 1 
cc, = gro (5 zs sin Pow) sin(a + 8), (A6) 
1 
CC = ST 3 sin? Ow sin(a + 8), | cca = m di sina + p cosa), (A7) 
Q11 = ccl cos? 0, + cc2 sin? 6, + 2cc3 sin 0p cos Op, (A8) 
Q1» = (cc2 — ccl) sin p cos 0; + cc3(cos? Oy — sin? 6p), (A9) 
Q»1 = (cc2 — ccl) sin p cos 0) + cc3(cos? 0j — sin? Op), (A10) 
Q22 = ccl sin? 0 + cc2 cos? 6, — 2cc3 sin 0p cos Oy, (A11) 
B; E Bi(p, mg,msi.)|yz =m? , (A12) 
B$ | = B;(p, m3, ms) ssa; Bi * = B;(p, m3, vg la. (A13) 
Cg! = Pai Mg» mz, mg), C mes Cij(—pi, —k, mg, Mg , m), (A14) 
C$! = Cij(— po, k, Tr» mg, mz), of! > Ci; (—pa, k, mg, m; , mi. ), (A15) 
C = Ci (— P2, —pn, mg, mi, Mg) Cs d = Ci (-pi pn, mg, mi, m; ), (A16) 
Di; H D(-— pi, —k, py, mg, mg ,, mi, mi), D} 17 = Dij(—pi Ph, pa, mg, m; ,, mj, mag), (A17) 
Di HU = Dij(—po, k, —pn, mg, mg, , mg, , mg), (A18) 
where B;, Cj; and Dj; are loop functions defined in [31]. 
After ge is factored out, the renormalization constant of b-quark can be expressed as 
167 
2 2 

6Z,=Cr) ABI, | 6Zg — Cp V, AFBI, (A19) 


I=1 I=1 
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where Cpr = 4/3. The contributions of the self-energy diagrams of b-quark propagator can be written as 


2 2 
DE =Cr) APB?" -6Z,, E$ =Cr)_ APB}! - 62Z,, 


I=1 T=: 
2 2 
De =Cr) | APB’ -ôZr, Y —CRM.ADPB|!-ÓZm. 
I=1 I=1 


C;, appeared in Eq.(2.5) are given by 


3 ama arya a a ^ 1 ^ 
Cs = 5 L7 58015 $05! 205, — 1/2 — mjOS AT [8 SCR AT [8 
I=1 
Bx ; shed ; 
-3 Ciz! Cs; + 205, — 1/2 — m505 DAT JE 30% AT f) 


ÔM, 1 


1 
Gtp 


R AI 1 
6Z,/8 —-óZg[t — X5,/8 — X JE (5641 + 50Zn + 


Mb 


2 
2mw cos B 4 e 2.4 ~ 3 
LL Qr;la4maC0 P Bp, /8 + 37408 P BEES 5s Do "BS 


Mp sin Q Fi 
EE 1 
e = $5605! + $2, + 22! — 1/2 - m20$ DAT [8 + C24 AT [8 
qd 
3 Wale Wale C c ; 1 n 
-3 Cis! + ICs) + 205, — 1/2 — m509 DAT [E 304 AT /t} 
i g 2 | 1 ó po 
jg E E A Lq m m n PS e 
2 2 my Ss t 
2 
2mw cos B 4 Pe Re an d IJ pR je, 2 hIJ pR 
tma 2a Og Bry /8+ Zm% Byj/t+ 5mgDo Bi 


2 

3 " " xd . I 

Cs = X L7 57405, + 1+ 2m508 AT [8 — (801) + 835. + 2054) AT /8) 
I=1 


4+26Z1,/8 + 203,/8+ (ZL + 6ZR 4 EUR \/8+ mind 


ó 2 : 8 
mb, ,. , 2mw cos 8 NS Qrz{—=mgCs  BE,J$ 
I 3 


J=1 


1 1 
*3ma( DG uu pur E Die Br -3m;(Df,7 — Di.'”) Bry + 3a Dis BIs}. 


2 
3 a a â 1 â 4 â 
Cs = 220-2408 +1+2m20¢ "AT /3 — 3 6Cts! + 8C% + 20%) AE /8} 
I=1 


OI 


2 a 
+28ZR/8 +25$/8 + (ZL +ôZR +2 dri Cong 
Mb 


2 

ae ; 
X Qr;i-3maC0 ™ Br, /8 
I,J=1 


)/8 + 


Mp Sin a 
1 1 
+3mg(Do My Dg — DLE zs 3m; (Dj - Dia) Bry ca zma Dis Br), 


2 
3 " z a T . " 
Cy = » (^75 0408 +1+2m2C§ P) A /E— 3 iC ta! T4085 2087 ATA 
I=1 


A ^ omy... 2mw cos B : 8 IJ T 
20Zg/t + 2X, /É + (ZL -6Z 2——)/t + ———— c mic gb 
T20Zn/ n/t * (Zr + óZn + "m )/t + TETTE 2 Qi 3'naCo 1j ft 


14 


(A20) 


(A21) 


(A22) 


(A23) 


(A24) 


(A25) 


(A26) 


1 1 
*amsDfs" Bry + 3m (Dja E Di3"”) Bry + zma(Do RE Di = D{,)BE}, (A27) 


2 
3 iy dus ; à 
Cio = X 4-350403 +1 + 2m306 Ar [E- S01! + EOS! + 2C54") AF /E} 
I=1 


^ n ómp, ,, 2mw cos 8 2 8 IJ n 
26Z, /i+ 20% /f + (6Z adu p er EE — m;Cf '! BE, jt 
+2821 [b + 29, [Ee (621 +8Zr e 27) fie =D) Qra{— gag " Bol 


1 
Ens D1, BR, + 3mq(Diy!! — DISP )Bfs + 


E mg(D§ + Di, — DS)BE. (A28) 


Since we have neglect the b-quark mass throughout this paper, C1,2,7,8,11,12 are irrelevant to our result and we do 
not present their explicit forms here. 
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